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ABSTRACT 

Hawking's radiation effect of Klein-Gordon scalar field, Dirac particles and Maxwell's electro- 
magnetic field in the non-stationary rotating de Sitter cosmological space-time is investigated by 
using a method of generalized tortoise co-ordinates transformation. The locations and the tem- 
peratures of the cosmological horizons of the non-stationary rotating de Sitter model are derived. 
It is found that the locations and the temperatures of the rotating cosmological model depend not 
only on the time but also on the angle. The stress-energy regularization techniques are applied 
to the two dimensional analog of the de Sitter metrics and the calculated stress-energy tensor 
contains the thermal radiation effect. 

Subject headings: Rotating non-stationary de Sitter cosmological model; Hawking radiation; generalized 
tortoise coordinate transformation 



1. Introduction 

The black hole thermal radiation technique of 
quantum f ield t heory was initially discovered by 
iHawkingl (I1974I . [19751 The quantum thermal ra- 
diation due to black hole has been investigated 
by different authors in diff erent types of space- 
time, such as Kerr- Newman ( Damour and Ruffini] 



1976 ). Kerr- Newman-Kasuya ([Ahmed and Mondal 



I995I ). Kerr dZhao. e^.a l."l983V Vaidya-Schwarzchild 
de Sitter (Dai. e^.a/Jll99 3) and NUT -Kerr-Newman - 



de Sit ter ( Ahmed .111 99 ll ) space-time. Zhao and Dai 
(|l992l ) proposed a suitable method known as 
generalized tortoise co-ordinate transformation 
(GTCT) which can determine the exact value of 
both the location and the surface gravity of the 
horizons of non-stationary black hole. The idea of 
GTCT is to reduce the scalar field equations, like 
Klein- Gordon or Dirac, in a black hole background 
to a standar d wave equation near the horizon. 
Li and Zhao I (|l993f ) investigated the Hawking's 



effect of a Dirac particles in a non-static case. The 
study of Hawking evaporation of Dirac particles in 
a non-stationary axisymmetric black hole is com- 
plicated because there is no unique method for the 
separation of variables for the radial and angular 
parts involved in the Chan drasekhar-Dirac equa- 
tions ( Chandrasekhar ] ll983l) of non-stationa ry ax- 
isymmetric space-time. IWu and Cai.l (|200ll ) stud- 
ied Hawking's effect on various non-stationary 
space-times, rotating Vaidya, non-rotating Vaidya 
de Sitter (2002a) and rotating Vaidya-Bonnor 
(2002b). For an evaporating black hole, the vac- 
uum expectation value of the renormalized en- 
ergy momentum tensor determines the location 



and temperature of the event horizon (jHiscock 



1981^ . It is well known that the de Sitter line 



element is a maximally symmetric curve space- 
time. The non-rotating de Sitter space is con- 
formally flat Cabcd — space-time with con- 
stant curvature Rgbcd = {A/^){gacgbd - gadgbc) 



(jHawking and Ellis .1119731 ) 
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The aim of this paper is to investigate the 
Hawking's radiation effect of Klein-Gordon scalar 
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field, Dirac particles and Maxwell's electromag- 
netic field in non-stationary rotating de Sitter cos- 
mological space-time. The non-stationary rotat- 
ing de Sitter solution is algebraic ally special i n 



Petrov classification of space-times ( Ibohal, 2009() . 



whereas the stationary rotating de Sitter is Petrov 
type D possessing a geodesic, shear-free, expand- 
ing as well as non-zero twist null vector which is 
one of the repeated pri ncipal null di rections of the 
Weyl curvature tensor ( IbohalJ[2005f) . For the sake 
of completeness of the presentation, we cite that 
the non-rotating stationary (Hawking and Ellis 



lis] 

m 



19731 ) and non-stationary de Sitter ()IbohalJl200 
models are conformally flat {Cabcd = 0) solutions 
of Einstein's field equations. It is worth to men- 
tion that the rotating de Sitter models, stationary 
as well as non-stationary, possess four roots of the 
polynomial A = r^ — A(M)r'*/3-|-a^ = 0, which de- 
scribes the singularity of the models, whereas most 
known black holes (rotating and non-rotating) 
have two roots of the polynomial. It is noted 
that the non-rotating de Sitter spaces, stationary 
and non-stationary, with a = have two roots 
of the polynomial A = 0. Gibbons and Hawking 
(1977) have suggested that the de Sitter cosmo- 
logical event horizon has formal similarities with 
a black hole event horizon. Therefore, it would 
be interesting to investigate the radiation spectra 
due to the Klein-Gordon scalar field, Dirac par- 
ticles and Maxwell's electromagnetic field in the 
rotating de Sitter space having four roots of the 
polynomial describing the cosmological horizons. 

The outlines of the paper are as follows: In Sec- 
tion 2, we reintroduce the non-stationary rotating 
de Sitter solution with cosmological function A(u) 
in Newman-Penrose (NP) formalism to be used in 
the study of scalar field equations in the next sec- 
tions. The NP spin coefficients of de Sitter solu- 
tion are helpful to study of the second order wave 
equations as seen in the remaining sections. In 
this section we also introduce an improved GTCT 
transformation in the non-stationary rotating de 
Sitter cosmological background having four roots 
of the polynomial, which describes the singulari- 
ties of the solutions. The locations of the horizons 
near r = rj,j = 1,2,3,4 satisfying the null sur- 
face condition are given for discussion in the next 
sections. The Klien-Gordon equation of scalar 
field is analyzed in Section 3. The adjustable 
parameters Kj appearing in the GTCT near the 



horizons are found for the Chandrasekhar-Dirac 
and Maxwell's equations in Section 4 and 5 re- 
spectively in the rotating de Sitter background. 
In Section 6, we recast each second order equa- 
tion to a standard wave equation of the Klein- 
Gordon scalar field, Dirac particles and Maxwell's 
electromagnetic field, and study the thermal ra- 
diation spectra of the fields by expressing the 
exact values of Hawking temperature. Section 
7 deals the calculation of the stress-energy ten- 
sor for two dimensional rotating de Sitter met- 
ric, and it is found to contain the thermal radi- 
ation in the time-time component guu of the met- 
ric tensor. Section 8 deals with the conclusion 
and remarks of the results obtained in the paper. 
The presentation of this paper is essentially based 
on the Newman-Penrose spin- coefficient formalism 
( Newman and Penrose] 1962 ) in (-1-,—,—,—) sig- 
nature. 

2. Non-stationary rotating de Sitter so- 
lution and generalized tortoise co- 
ordinate transformation 

Here it is convenient to r eintroduce th e non- 
stationary rotating de Sitter (|lbohalJl2009l) in NP 
formalism for the study of Hawking's radiation ef- 
fect of Klein-Gordon scalar field, Dirac particles 
and Maxwell's electromagnetic field. The line ele- 
ment, describing non-stationary rotating de Sitter 
solution with a cosmological function A(u) in the 
null co-ordinates system (u, r, 9, 0) is 



1 - 



3^2 

2aK{u)r^ 



3i?2 

-R^de^ - {{r 

-Aa^ sin^e} R 



du^ + 2du dr 
'^9 dudcj) - 2a sm^9 dr < 

2\2 



sm 

2 



where A = - A(u)r^/3 + 



d(t)^ 
R 



(1) 



la cos ( 



^cos^6' and u is retarded time co- 



ordinate. This line element will represent rotating 
stationar y de Sitter so lution when A(u) becomes 
constant (jlbohal.1 120051 ). Also if one sets the rota- 
tional parameter a to zero, it will describe t he non- 
stationary non-rotating de Sitter solution ([Ibohal 
2009). Here it is convenient to express the null 
tetrad vectors explicitly for further use as 



a sm 



A 

2R^ 



2R'- 



■asm'^9Sl, 



2 



ma = - 



(2) 



Here £a, ria are real null vectors and is com- 
plex with the normalization conditions tan"" = 1 = 
—ma'm"' and other inner products vanish. The di- 
rectional derivatives are defined by 



1 



-I- 



V2R sm0 ^' 

Then, the NP spin-coefficients for the metric (1) 
are given as follows: 

K = (T = A = e = 0, 
_ 1 _ A 



a = TT — /?, /3 = 
ia sm 6 



cote 



V2RR 



2V2R\ 
ia s'ml 

T = — - 



(4) 



2A(w)r2 



)i?-A}, 



6V2RR'^ 



ia sin 9 A(u),„ r** 



From these NP spin-coefficients, it is to mention 
that the rotating non-stationary de Sitter model 
possesses null geodesic affinely parameterized (k = 
e = 0), shear free [a — 0), expanding {9 ^ 0) and 
rotating (w^ ^ 0) null vector £a- 



2 1, .n a'^cos'^9 



(5) 



By the virtue of Einstein's field equations Rat — 
(1/2) R gab = —KTab, the energy-momentum ten- 
sor describing matter field for the space-time (1) 
is given as 

Tab = H*£aib + '2p*£(^anb)+2pm(^arrib) 

+2uj£(^amb) +2uj£(^a'mb), (6) 

where the quantities are 

^4 



P = 



GKR^R' 



r{2rA(w),„ 4-a^sin^6lA(w),„„ }, 



KR^R^ 



A(w), 



KR^R- 



6V2KR^R^ 



-{r^ + 2a^coa^9}A{u), 
{R~3R)Aiu),u, 



(7) 



with the universal constant K = SttG/c^, /i* de- 
scribes null density arisen from the non-stationary 
property of the solution with cosmological func- 
tion A(u). p* and p are the density and pres- 
sure of the matter distribution which does not 
involve derivatives of cosmological function A(u). 
Here, u> represents the rotational density deter- 
mined by the rotational parameter a coupling with 
the derivative of A(u). However, uj will vanish 
when either a is zero for non-rotating, or A(u) be- 
comes constant for stationary models. The energy- 
momentum tensor (6) satisfies the cons ervation 
equat ions, expressed in NP formalism (jlbohal 
20091) 



T'^ = 0, (8) 

which shows that the non-stationary de Sitter met- 
ric (1) is an exact solution of Einstein's field equa- 
tions. 

We shall discuss the nature of time-like unit 
vector — {\/2)^^{£a + n°')- By knowing the na- 
ture of time-like vector field of the observer, we 
can explain the physical properties of a matter 
distribution of rotating line element i.e., whether 
the field is expanding (Q = u"-^ ^ 0), shearing 
(fafc 7^ 0), rotating {ujab 7^ 0) or accelerating 
{ii = Ua-fiU^). The expansion scalar Q and acceler- 
ation vector iia — Ua-.bU^ are obtained as follows 



e 



V2i?- 



2A{uy 



f 2A(u)H A -1 



2i?2 

a2sin^6' 

~ 2V2R^R^ 



{Rma + Rfa"-^ ■ (9) 



where Va — -^i^a — "n-a) and A = — A(u)r^/3 -f 
a^. Then the shear Oab and vorticity ujab tensors 
are found as 



CTafa = U(a;f,) " U(aUb) " ^©(Safc " UaUb) 

-4(r-^)" 



1 r3rA 



3\/2i?2 



i?2 



: {va 



X \VaVb - m^amb) 



1 



V2RR^ 



X 2aisin0(2r + iacos( 



3 



zsin 9A{u) 



6 



V(amb) + {c.c.} 



U[a:b] - U[aUh] 

ia sin C 



2ai cos 9 



A(u),„r4 



6 



W[a"^&] 



, . m cos f / 



— W[amb],(10) 



which are orthogonal to i.e., aabu^ = and 
oJabU^ = 0. Here {c.c.} denotes the complex con- 
jugate of the previous term. The expansion scalar 
Q and acceleration iia in (9) of the solution de- 
pend on the function A(w). The presence of rota- 
tional parameter a in the rotating metric (1) shows 
directly the existence of vorticity tensor in (10). 
When the rotational parameter a sets to be zero, 
the vorticity tensor ujab will vanish, indicating the 
nature of the non-rotating solution. 

The non-stationary rotating de Sitter solution 
(1) has singularity when 



A = r2-r''A(?i)/3 + a2 



0. 



(11) 



This polynomial equation has four roots r+ -|- , 
r+ - , "T- + and r They are found as 



-±(±) = ±\/ ± V9 + I2a^m }. (12) 



Now let us denote these roots r+ + , r+ _ , r_ 
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r as ri, r2, r^, r4 respectively for simplicity, 

Tj, j = 1,2,3,4. They satisfy the following rela- 
tion 



(r - ri)(r - r2)(r - r3)(r - r4) 
3 r o r*A{u) 



A{u) 



«1- 



(13) 



Then each root will describe the location of the 
cosmological horizons for the observer and asso- 
ciates an area of the horizon at each point at 
r = r,, j = 1,2,3,4, 



\/geeg4,4, dO dcj) 



JO 

2 I „2l 



The horizons are necessarily null surfaces r 
r{u, 9) that satisfy the null surface conditions 



The generalized tortoise co-ordinate transforma- 
tion (GTCT) technique gives the location and 
tempe rature of hor i zons of non-stationary black 
hole (IWu and Cai.l I2OOII) . Using GTCT, one 



can reduce Klein-Gordon or Chandrasekhar-Dirac 
equations to a standard form of wave equa- 
tions near the horizon by introducing common 
tortoise co-ordinates transformation by = 
r + {2K)^^ln{r — rn) in a non-stationary or sta- 
tionary space-time, where k is the adjustable pa- 
rameter and th is the event horizon. We assume 
that the locations of the de Sitter horizons Vj are 
functions of retarded time co-ordinates u = t — 
and angle 9. 

Here, since the rotating de Sitter cosmological 
space has the four roots, we would like to introduce 
an improved GTCT as follows: 



1 



2nj{uo,do) 

u — Uq, Of, = 



ln[r — rj{u, 0)], 



(14) 



where Vj = rj(u,6) are the locations of the hori- 
zons and the adjustable parameters to be 
determined. All parameters Kj , uq and are con- 
stant under the GTCT. From the above equation 
(14), we have the following operators 



d _ d 
dr dr^ 
d _ d 

d9~d9^ 

d _ d 
du duit 



1 



d 



2k. j (r — rj ) 9r, 
r, d 



2Kj{r — rj) dr^ 



(15) 



d 



2Kj{r — rj) dr^, 



where fj — drj/du, Vj — drj/dO. Now applying 
the GTCT in (14) to the null surface condition 
gabQ^pQ^p _ Q taking limit as r — >• rj{uo, 9q), 
u ^ Uq and ^ 00, where Uq and 0q are the initial 
states of the cosmological horizons, then we get as 



+ 2(r| + a^)rj + a^s\T?0ofj 



(16) 



The above equation holds only when the content 
of the square bracket is equal to zero giving the 
location of horizons of rotating de Sitter cosmo- 
logical model as 



g'^'daFdbF = 0. 



+ 2{r] + a^)rj + a^sm^0of] + r] = 0, (17) 



4 



where Aj = r'j—A{uo)rj/3+a'^. Here fj = drj/du 
is the rate of horizons with respect to time and 
fj = drj/dO is the rate of change with angle. 
Equation (17) shows that the shape of cosmologi- 
cal horizons depend on time and angle due to the 
GTCT transformation. Equation (17) has four 
roots as in (12) and they are given below: 



^{3 + 6f,±VTI^}],(18) 



where = 9 + 12a^K{u) + iH and H = I2fj{l + 
Tj) + 4A(M){r| + 2a^fj + a^f^ sin^ 9}. These roots 
are due to the GTCT in (14) and will reduce to 
those of the non-stationary rotating de Sitter cos- 
mological model given in (12) when rj — fj — Q in 
(18). 

3. Klein-Gordon equation 

The wave equation describing the dynamical 
behavior of scalar particle of mass /iq is given by 

(5aV^5"'4$) - Mo'f = 0. (19) 



For the de Sitter line element (1), this equation 
reduce to 



2a- 



92 



-2a 



32 

H 

'r 96*2 sm 

d , d 



dudr dud4> 
1 92 



2r 



AK{u)r^\ d 
dr 



$ = 0. 



(20) 



Under the GTCT and taking the limit r — > 
rj(uo,Oo), as u — > uq, — > Oq, the above equa- 
tion can be written as 



-{• 



2A(uo)r3 



3 
92$ 



2f 



2A,- 



u -f 92$ 

+2(r2 + a2) — - 2fj—— 

-2«(i + ^.)iSr-2{(^,' + «') 



92$ 



- {2rjfj + a^sin^Oofj 
0. (21) 



In deriving the above equation, we use L' Hospital 
rule to deal with infinite form of 0/0 and the value 
of the limit is given as 



lim 



A + 2(r2 + a'^)fj + a'^sm^Oaf'^ + f'^ 



2A(iio)rf 



'^rjfj 



(22) 



The standard form of wave equation near the hori- 
zons can be obtained by adjusting the parameter 
Kj for each rj in (21). The wave equation (21) 
shall be discussed later in Section 6. 

4. Chandrasekhar-Dirac equations 

The four coupled Chandrasekha r-Dirac equa- 

tions in Newman- Penrose formalism (Newman and Penrose 



19621 ) are given by 



p)Fi- 


f(^ + 


TT — 


a)F2 


= ifJ-oGi, 


l)F2 


+ (<5h 




t)Fi 


= ipaG2, 


P)G2 




- 7f — 


a)Gi 


= iMo^2, 


l)Gi 


-(^ 


f ^- 


-f)G2 


= ip-oFi 



[D + e 



(23) 

where no denotes the mass of the Dirac particles 
and Fi, F2, Gi and G2 are the four components 
of the wave function. By substituting the spin- 
coefficients from (4) in (23), we have the following 
equations 

{dr + i)Fl + ^^"^2 = iMoGl, 

A „ I „ 1 , I iasin0,„ „ 

kD+F2 + -^(L+ + = — )Fi = iuoGa, 

2i?2 ^i?^ R ^ ^ ^" ^' 

-^L+Gi + (dr + ^)G2 = ipoF2, 



la sm f 
R 



■)G2 = i^aFi, 



(24) 



where operators are given by 

= dg + ia sin 6* 9„ 4 
= da — ia sm6 du — 



smf 
i 



cote 
cot 9 



D+ = 9,. - A-i|2(r2 + a2)a„ 



2a ds 
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(25) 



We define Fi, F2, Gi and G2 by 



TP n ^1 

i^2 = /2, G2 = ||. (26) 



Then equation (24) becomes 
dfi 



dr 



L f2 = iiioRgi, 



dr 

AD+gi + L-g2 



-2ifxoRfi. (27) 



We may consider with a pair of components /i , /2 
or gi, g2 only. In order to investigate the Hawk- 
ing's radiation, wc will consider equation (27) near 
the horizons only. Under the transformation (15) 
the limiting form of equation (27) near the hori- 
zons are obtained as follows 



dfi ., . . . ^.<Jj2 
__(,._, a,. sin0o)^=O, 

|rj + fjia sin^oj^^ + |Aj 

+2'^.(r-| + a^)}|^=0, (28) 



9/2 



and 



{A,+2.,(.|+a^)}|i^{.', 



—rjia smUo >— — 



= 0, 



(29) 



after taking the limit as r ^ rj{uo,6o), u ^ uq 
and 6 ^ 6o- If the four derivatives 
^ and of the above equations (28), (29) arc 
nonzero, then the condition for non-trivial solu- 
tions for /i, /2, gi and g2 is that the determinant 
of its coefficient will be zero. This gives the equa- 
tion: 

Aj + 2(r| + a^)fj + a'^sm^Oof] + r| = 0, 



where Aj = ^ A(?io)r^/3 + . This equation 
is same as the equation for the null surface (17) 
given above. The equation (28) can eliminate the 
crossing term of first order derivatives in second 
order Chandrasekhar-Dirac equations near r = rj. 

A direct calculation gives the second order 
Chandrasekhar-Dirac equations as follows 

AD+^ + L-L+h - ^tilR^h - *A^oAffi 



-a Ho sin + 



ias\n.Or^k{u) ( df2 



d_ 

dr 



(AD+/2) + L+L-/2 - 2/X^i?2/2 

= -a/io sin 6gi - Z/UoS'2- (30) 



Using the GTCT in (30) and after some tedious 
calculations, the limiting form of above equation 
when r tends to rj{uo, Oq), u approaches to uq and 
6 goes to 60 leads to 

{2A,+2(rf + a2)r, + '"^'+^'"^'''^' 



3/^ 



-2r 



d(l) 



+a sm 



_ 2A(Mo)rj 
dr^du^ 3 
+rj + Vj cot 9o — ia cos 6ofj 

2 2 msin6'oA(wo)?'^ 
+a sin 9orj + ^rjfj 

fj+iasm9ofj ^ dfi _ 
A,+2f,(r2+a2)/ar. 



(31) 



and 



r,- -I- 2rjfj 



{2A, + 2(r|+a^)f, + I^ 



-2r 



2A(Ko)r| 

^^2 

+ rj +fjcot6o 



r _ 2A(wo)r| 



+ 



6 



0orj +aW0or,}^ =0, (32) 



~ia cos y^rj 

here the first-order derivative term dr, /2 in equa- 
tion (31) is replaced by 9r,/i- In order to reduce 
to a standard form of wave equation near the hori- 
zons, we adjust the parameters Kj such that they 
satisfy 



-{' 



2A(uo)? 



+2Aj = r| + + a^sin^^orj, (33) 



and using (17) we obtain the surface gravities of 
the horizons r — as 



2A(uo)r^ 



Q(l + 2f,) 



2.7 



(34) 



where Q — 



a'^sin^Oorj with j = 1, 2, 3, 4. 



Here we observe the involvement of the de Sitter 
cosmological function A(m) at m = uq in the ex- 
pression of surface gravities. These surface grav- 
ities Kj are the parameters for determining the 
temperature of the thermal radiation spectrum for 
scalar fields. 

5. Maxwell's equations 

In this section we shall analyze the Maxwell's 
equations in non-stationary rotating de Sitter 
space-time. The four coupled Maxwell's equations 
for electr omagnetic field in Newman-P enrose for- 



malism (jNewman and Penrose 

by 



19621) are given 





- S(j)o 


6(t)2 - 


-V01 




- V0O 


D<j)2 


-Hi 



(35) 



(tt - 2a)0o + 2p(j)i 
-K(j)2 + 27roJi, 

+(t- 2/3)02 + 2^0^2, 
in - 27)00 + 2t0i 
-0-02 + 27ro J3, 
— A0O -I- 27r0i 

+ (p-2£)02 + 27roJ4, 



where ttq = (22/7) and D = i^da, V = n^da, 
and 6 = ra°-da are the directional derivatives and 
Ja is the tetrad components of the current vec- 
tor given by Ja = jfiZj^. In the above Maxwell's 
equations, by inserting the spin coefficients (4) 
and substituting 0o = ^0, 0i = {V^R)^^Fi and 



02 = 2 ^F2 to the above equation (35), we 
may write them as follows: 



dr RJ ^ 
I aisinf 



^0 



R 

ai sinf 
R 



Lo 
Fi + 
F2 + 



R 

ai sin i 

a(d+- 



Fo = 0, 
Fi=0, 
i)Fo = 0, 

rJ 

(36) 



where operators are defined by 



Ln 

Dt 



de — ia sin 6 du 
de + ia sin 9 du 

a^-2A-i|(r~ 
2A(u 



— nr 1 — 



sin 9 

i 

sin 9 

a')a. 

-)}■ 



n cot ( 



d^ + n cot 9, 
i + add, 



(37) 



Here n is positive integer. The first order 
Maxwell's equations near the horizons are as fol- 
lows 



rj — rjta sill 9o 



dFp 

dr.. 



dFi 
dr^ 
\ dFo 



= 0, 



+iasm9orj j-^— = 



0, 



(38) 



and 



rj — rjiasmOo 



dFi dF2 

dr^, dr^, 
\dF\_ 



= 0, 



1 dF2 

-t-iasin^o'^i (-^ — = 0. 
J or. 



(39) 

We assume that the derivatives -^Fq, ^-^1 and 
■^-F2 are non-zero; and similar to the Dirac equa- 
tions (28, 29), the condition for non-trivial solu- 
tions for Fq, Fi and F2 are that the determinant 
of the coefficients of each (38), (39) will be zero. 
That provides the horizon equation like the Dirac 
equations or the null surface condition (17). Then, 
equation (39) can eliminate the crossing term in 
the second order Maxwell's equations near the 
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horizons. The second order Maxwell's equations 
are given by 

i+i,F„ + -(AD+F„)-2«^=0, 



^2ahm^9Fi _ 4A(u)r3Fi _^ 2AFi 



RR 



3R 



RR ' 



^ + LoL, F, + — 

+2rM = ai^^^(^MuyF, 

du 3R ' ^ ' 

The second order Maxwell's equations become 



[-{• 



2A(wo)d 



2r 

9r2 ^ dr^dO^ 



2rj{l+rj) 



fj cot 60 + 2ai cos Oo'Tj 



dFp 
5r. 



(41) 



\-Kn I' 



2A(uo)r: 



}+2A,- 



+2fj{a'^ + r^) 



|-^ + 2r,r 



d F 

-2a(l + fj)—^ - 2{r'] + 



+a^sin^^o'^j) 



d^Fi 



+rj + 2rjfj + fj cot Oq 



r 1 r 2A(uo)r:^ 1 
[-{r, ^+2f,r,}+2A, 



dFi 
9r* 



= 0, (42) 



+2f,(a2 + r2) 



9r2 



2r, 



d^F2 



2r,- 



4A(Mo)r| 



+Arjfj + + cot + a^sin^^o^j 



+2(rj + ia cos Oo)'fj 
fj + ia sin ^o^^i 



Aj + 2{r^ + a^)fj 



ai sill f^nA((/|))r| 
3 



(43) 



after taking the limit as r — >■ rj{uo,9o), u — )• 
Wo and 9 ^ 9q. By reducing the second or- 
der Maxwell's equations under the transformation 
GTCT, the parameters Kj near the horizons are 
obtained as follows 



2A(«o)rf 
3 



2rjrj 



g(l + 2f,) + 2r| 



(44) 



which are the same as (34) above - the sur- 
face gravities obtained from Chandrasekhar-Dirac 
equations. 

6. Hawking temperature and thermal ra- 
diation spectra 

The Klein-Gordon equation (21). Dirac equa- 
tions (31), (32) and the second order Maxwell's 
equations (41), (42) and (43) could be recast into 
a combined form near the horizons Vj as follows: 



52^- ^ r 92^ 
^-2|ifi 



^2* 



924- 



dr*d9* ' dr^fdcf) dr^.du^ 

(i^2 + ii^3)|^}=0, (45) 



where 



(46) 



where Q = rj + + a2sin^^o fj ■ The equation 
(45) is the standard form of wave equation in the 
non-stationary rotating de Sitter space. This wave 
equation includes the Klein-Gordon, Dirac and 
Maxwell's electromagnetic field equations having 
different conditions of the constant terms. 

For example, when ^ = ^ for Klein- Gordon 
scalar field, the equation (45) yields the following 
constants as 



{2rjfj + a sin ^o^^: 



2Q 

+ cot 6ofj + Tj 



8 



0. 



(47) 



For Dirac particles when 4' = /i, the constants 
are given by 



1 r 
2Q 



— cot O^r-j — Vj 



- a^sin^^or, 

3 



3{A,+2f,(r2+a2)} 



K3 = — 
^ 2Q 



cos 6'f 



rjr|sin0oA(^io) 



^^3{A,+2f,(r2 + a2)} 



(48) 



and for 5" = /2, we get 



where ifi, VI K2 and i^a are considered as real 
constants under the generalized coordinate trans- 
formation (14). One can view that equation (45) 
is a second order partial differential equation rep- 
resenting wave equations near the horizons since 
all the coefficients Ki, ^j, K2 and K3 are taken 
as real constants when r — >■ Vj {uq, Oq)^ w — > uq and 

Now following Wu and CaiJ ( 2001 ). we may 
separate the variables in (45) for the analysis of 
the field equations as 



* = R{r*)Q{e*)e 



lymcp—uju^ 



(53) 



The radial and angular parts of equation (45) are 
respectively given by 



1 

K2 = — 
2Q 



2A(uo)rf 



-•fj cot 9q — a^TjSm^Oo 



K3 = —acosOorj. 



(49) 



Similarly, for Maxwell's electromagnetic field = 
Fq, the equation (45) provides the following con- 
stants: 



K2 = 

K3 = 
for * = Fi, 



1 



r,(l + r,) 



3 



r 



,2/ 



J r.j cot t^o a svcCOi^rj 

^ \ 2 2 

a cos O^fj 



1 



(50) 



K3 = 
and for * = i^2, 



K2 = -^{a^sin^eofj 

+2rjrj + Tj cot^o), 

0; 



(51) 



K2 



1 

Q 



Tj cot Oq 



2A(uo)^ 



a^sin^^ofj 



2 . 2 

a'^ sin^ 6*0 r j A (uq ) 



6{A, + 2(r2+a2)r,} 



Ka 



cos t^o^'j 



Tj sin ^?oA(uo)r 



6{Aj + 2(r 



(52) 



8 =Ae 
52 i? 



2{ir7ir2j 



1K3 - Xo} 



Or, 



-0, 
(54) 



where A is a real constant of the variable separa- 
tion, Kf) ~ K2 — KiX. The ingoing and outgoing 
wave solutions to (45) near the horizons are given 

by 



^out 



{r- 



X (r 



r > r 



' 



(55) 



where ^'JJ^ represents an incoming wave and is 



analytic on r 



denotes an outgoing 



wave having logarithmic singularity on the hori- 
zons, but not analy t ic at r = r,. Fo l lowin g 
Damour and Ruffini J (|l976l ) and ISannanI (|l988t ). 



it can analytically continue from the outside of the 
horizon r, into its inside as 



(r - rj) 



(56) 



to 



X {t 



0" 



■(57) 



when (r < rj). The relative scattering probabili- 
ties at r = Tj are given by 



- {niilj —uj — K^ ) 



(58) 
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where m is the azimuthal quantum number, w is 
the energy of the radiating particles and flj is the 
angular velo cities near r = r^. Accordin g to the 
suggestion of iDamour and RuffiniJ (|l976l ) and ex- 
tended by ISannan ( 1988 ). the thermal radiation 
spectrum is given by 



r /'I 
(7V^) = {cxp(- 



iQj — cj — K3 



with 



27rl 



2A(Mo)r| 



0(1 + 2rj) 



2r] 



(59) 



(60) 



which is the Hawking's temperature defined by 
Tj = Kjllix. This determines the temperature 
of the thermal radiation near the horizons Vj 
due to Klein-Gordon scalar field, Dirac particles 
and Maxwell's electromagnetic field in the non- 
stationary rotating de Sitter cosmological space- 
time. Since the surface gravities Kj are same for 
these field equations, the temperatures of thermal 
radiation spectra are same as in (60) depending on 
the de Sitter cosmological function A(-uo) near the 
horizons. However, it is observed that the thermal 
radiation spectra are different for different wave 
equations as the constant Kj, are not the same for 
all the scalar fields discussed here. 

7. Description of the motion of rotating de 
Sitter metric in two dimensional space 

The line element of rotating de Sitter solution 
in Boyer-Lindquist co-ordinate can be written as 



2asm^0 



i?2 

{o^sin^eA - (r^ 



R^de 
) - Ajdtdcp 
a'f}d(b^ 



(61) 



where A = - A{u)r'^/S + and R'^ ^ + 
a^cos^d. In order to describe the space-time ge- 
ometry of the rotating de Sitter metric in two di- 
mensional space, we take — Ojir and d9 = d(j) = 
in the equation (61) and the two dimensional line 
element can be written as 



A 



R' 



(62) 



where R^ 



a? and / — AR 2. On the sym- 



metric axis, horizon is defined by gu = 0, whereas 



in general it defines the ergosurface. The two di- 
mensional line element could not describe the an- 
gular momentum and the transverse degree of free- 
dom are simply redshifted away relatively to the 
ones in the rt-plane. To calculate the stress en- 
ergy tensor for the metric (62), we transform to 
the double null form as 



f- 



dudv. 



(63) 



where u ~ t — r* ,v ~ t + r* , where r* is de- 
fined as tortoise coordinate given by dr* jdr = 
R^ / A. We obtained stre ss energy components 



(|Birrell and Davies.1 Il986l ) for two dimensional 



non-stationary rotating de Sitter space as follows 



T — T — — - 



1 



1927r 
4{AA,rr + 



A,l -4A{1 - 2K{uy} 



A2) 



R^R^ 
12A2r2 



T — 

11.1! 



1 

48^ 



A{l-A(w)r2} A(2rA,^ 



-A) 



i?2i?2 



R^R^R^ 



4r2A2 



R^R^R-^R^ 



T — — - 



1 



A,2-4A{l-2A(M)r2} 



1927r 
4(AA 



i?2i?2 



A2) 



12A 



2^2 



R^R^R^ 



R^R^R^R^ 



Tu = - 



1 

96^ 

8A(2rA,^-fA) 28r2A2 



A,2 -8 A(l-2 A(u) 

i?2i?2 



Ttr — 



R^R^R^ 

Trt = 0. 



R^R^R^R^ 



(64) 



From above, we observed that there is no constant 
terms involved in the expression r„„. The Hawk- 
ing flux can be obtained by considering the value 
of Tuu at the energy flux through the cos- 



mological horizon, where r-|_+ 
Since A = 0, Tu„ is given by 



3+-^9-12a2A(tt) 
2A ■ 



1 



1 



2A(«)r^ 



vv r— r4 



487r 



{rl+ + a2)2 



(65) 



which is exactly the negative of the Hawking flux. 
The temperature defined by Tj = Kj/2TT of the 
non-stationary rotating de Sitter solution at the 
cosmological horizons is given by 



T = - 



A{u)r 



TT l(r^+ + a2)2/- ^00^ 
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In the above equation, the negative sign denotes 
the radiation of the energy-momentum tensor flux 
is radiated at infinity. 

8. Conclusions 

In this paper we have studied the Hawking's 
thermal radiation effects of Klein-Gordon, Dirac 
and Maxwell's electromagnetic fields in the non- 
stationary rotating de Sitter background space un- 
der the GTCT transformation. It is interesting to 
emphasis some of the physical properties of the 
de Sitter solution (1). The non-stationary rotat- 
ing de Sitter model possesses the Weyl scalars 
-02 = = ^^4 l|lbohal.|[2009l ). These non- 



zero Weyl scalars determine the gravitational field 
of the space-time that the de Sitter solution (1) is 
an algebraically specia l type II \C nhr[d(-h](-^ (-'^ = 0, 



with ^/jQ = -01 = ( Chandrasekhar J 1983 ^] in 
the Petrov classification of space-time with a null 
vector £a. The nuU vector is geodesic, shear 
free, expanding [6 ^ 0) as well as non-zero twist 
(w 7^ 0) given in (5). Here the function A(m), de- 
scribing the non-stationary status of the solution, 
does not involve in the expressions of expansion 
9 and twist Cj in (5). This means that the phys- 
ical properties of this null vector 1°" are same for 
both stationary (with constant A) as well as non- 
stationary (with function A(u)) rotating de Sitter 
models, though they have different gravitational 
fields possessing different energy momentum ten- 
sors. 

The equations (17) and (34) above deter- 
mine the locations and surface gravities of non- 
stationary rotating de Sitter cosmological model, 
and their values depend on both the retarded time 
u and angle 6. Under the GTCT transformation, 
the second order equations of the Klein-Gordon, 
Dirac and Maxwell's electromagnetic fields are 
transformed to the standard wave equation near 
r = rj, and separation of variables can be done as 
in (53). It is observed that the constant K3 is ab- 
sent in the spectrum of Klein-Gordon scalar parti- 
cles in (47). Also K3 of Dirac and Maxwell's equa- 
tions will be vanished when the non-stationary 
de Sitter solution (1) becomes stationary with 
Tj = rj ^ or it has a zero angular velocity 
with a — 0. We have also seen the effect of 
the cosmological function A(u) in the tempera- 
tures of thermal radiation spectra for all the three 



wave functions under the GTCT transformation in 
non-stationary de Sitter background. For a non- 
stationary non-rotating line element with a — 0, 
the constants Qj, in both wave equations be- 
come zero, and the other constants will not vanish. 
For stationary non-rotating metric, the constant 
K3 will be zero for wave equations. The surface 
gravities for Klein-Gordon scalar field, Dirac parti- 
cles and Maxwell's electromagnetic field equations 
under the transformation of GTCT are the same. 

It is interesting to mention that the angular ve- 
locities flj given in (46) are due to the GTCT 
transformation in the non-stationary rotating de 
Sitter cosmological space: 



a(l + fj) 



a^sin^^i 



(67) 



The appearance of rj in this equation indicates the 
effect of GTCT in the rotating de Sitter space. 
The rotating de Sitter cosmological space-time, 
stationary or non-stationary, will have the same 
angular velocities ftj under the GTCT transfor- 
mation. When fj — 0, ftj will be those of rotating 
de Sitter space before the transformation as 



n = 



which depends on the roots r ~ Vj, j = 1,2,3,4 
of the polynomial A = 0. It is observed that the 
angular velocity for every rotating space-time, sta- 
tionary or non-stationary, has the same structure 
as in (68). But it depends on the roots of the 
polynomial. 

For example, the rotating Kerr black hole has 
the angular velocity ft — a{r'^ + a'^)^^ , where 
r — M ± y/{M'^ — a'^); similarly Kerr- Newman 
having r — M ±^y (M^ — a? — e^) with the charge 
e; Kerr-Newman-Vaidya (jlbohal and Kap il. 201ol) 
having r = M + f{u)±y/{{M + f{u))^ - - e^}, 
where f{u) is the Vaidya mass; rotating Vaidya- 
Bonnor, r = f(u) ± \/{P{u) — — e^(u)} with 
the charge e{u) and so on. From the study of 
Hawking's radiation effect of Klein-Gordon scalar 
field, Dirac particles and Maxwell's electromag- 
netic field in the non-stationary rotating de Sit- 
ter cosmological space-time under the GTCT, it 
is also found that the equation (45) could be con- 
sidered as a general standard wave equation to 
be tackle for every non-stationary rotating space. 
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The example of this situation is that three differ- 
ent wave equations, namely Klein-Gordon, Dirac 
and Maxwell could be recast into this standard 
form (45) with the variation in the values of the 
constants K2, K^. Also all these three wave equa- 
tions have the same thermal radiation tempera- 
ture (60) with the parameters Kj (33) or (44) for 
the non-stationary rotating de Sitter cosmological 
space-time. 

The difficulty in the separation of variables of 
the standard wave equation (45) may be seen that 
if one assumes 

* = i?(r*)e(r)*(</>)exp(-iwu,), (69) 

in (53), then the relative scattering probabilities 
at r = rj are given by 

> r,) 2 



^°^\r < r,) 



-if 3 



(70) 



where f3 is the separation constant Z = a-\-ifi. The 
equation (70) does not involve the angular veloc- 
ities VI j given in (46) or (67), which shows the 
difference from the scattering probabilities (58) 
with the separation of variables (53). The equa- 
tions (53) and (69) show at least two possibilities 
for the separation of variables in non-stationary 
space-times. However, it seems that the separa- 
tion in (53) may be a convenient one, because it 
shows the direct involvement of the angular veloci- 
ties VLj without disturbing the rotational character 
of the background space-time. 

The investigation of Hawking's radiation ef- 
fect of Klein-Gordon scalar field, Dirac parti- 
cles and Maxwell's electromagnetic field in non- 
stationary rotating de Sitter space with cosmo- 
logical function A(w) is not been seen discussed 
before. The work discussed here in the non- 
stationary rotating de Sitter model includes the re- 
sults of the following models: (i) stationar y rotat 



ingde Sitter when A(u) becomes constant (jlbohal 



20051 ). (ii) no n-stationary n on-rotating de Sitter 
when a = ( Ibohal.1 l2009l ) and (iii) stationary 
non-rotating when A(n) = constant and a — Q 
(jHawking and Ellis.|[l973l) . 
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